Dense Stellar Matter with Strange Quark Matter Driven by Kaon Condensation 
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The core of neutron-star matter is supposed to be at a much higher density than the normal 
nuclear-matter density, for which various possibilities have been suggested, such as, for example, 
meson or hyperon condensation and/or deconfined quark or color-superconducting matter. In this 
work, we explore the implication on hadron physics of a dense compact object that has three 
"phases": nuclear matter at the outer layer, kaon condensed nuclear matter in the middle, and 
strange quark matter at the core. Using a drastically simplified but not unreasonable model, we 
develop the scenario where the different phases are smoothly connected with the kaon condensed 
matter playing a role of a "doorway" to a quark core, the equation of state of which with parameters 
restricted within the range allowed by nature could be made compatible with the mass vs. radius 
constraint given by the 1.97-solar-mass object PSR J1614-2230 recently observed. 

PACS numbers: 97.60.Jd, 26.60.Kp, 26.60.-c, 06.30.Dr 



I. INTRODUCTION 



The recent observation of a 1.97-solar-mass (M@) neu- 
tron star, PSR J1614-2230, using Shapiro delay with 
Green Bank Telescope [l| , raises a highly pertinent is- 
sue on whether the non-nuclear degrees of freedom such 
as meson condensations (pions, kaons, and/or hyperons), 
quark matter (normal or color-superconducting) , and the 
multitude of other phases are relevant for the physics of 
stable compact stars. This observation, if unequivocally 
confirmed, will have far-reaching ramifications on sort- 
ing out the large variety of models proposed in the lit- 
erature for the maximum (minimum) mass of neutron 
stars (black holes). A particularly poignant example is 
the scenario of Bethe and Brown [2] [3], where the onset 
of kaon condensation inside the neutron star matter at a 
density p ~ 3po ~ where po is the nuclear matter density 
- keeps the maximum mass less than 2Mq, which seems 
to be consistent with the observations of well-measured 
neutron-star masses ~ 1.5M 0], beyond which, it is ar- 
gued, the stars collapse to black holes, giving rise to an 
enhanced number of light-mass black holes in the uni- 
verse. At a much more drastic level, there is also the 
issue of possible non- Newtonian gravity intervening with 
"soft" equations of state that may otherwise be ruled 
out by the observation [5|. This circle of issues was re- 
cently raised and analyzed, among others, by Lattimer 
and Prakash 161. 
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Unlike hadronic matter at high temperature for which 
lattice QCD provides model-independent control of ef- 
fective theories, cold baryonic matter at high density has 
remained a largely unchartered domain, more or less in- 
accessible by model-independent approaches. Because of 
the sign problem, lattice QCD cannot yet probe the rele- 
vant density and furthermore there are no known model- 
independent tools available for treating baryonic matter 
at high density. The abundant and accurate experimen- 
tal data on finite nuclei allow one to determine the equa- 
tion of state (EoS) for both symmetric and asymmet- 
ric systems up to the normal nuclear-matter density po- 
However, there is neither a reliable theoretical tool nor 
unambiguous experimental information to enable one to 
extrapolate the EoS beyond po in a controlled way. This 
conundrum accounts for the widely diverging model pre- 
dictions available in the literature. 

In this paper, we would like to address the EoS of 
compact-star matter, with the observation of the PSR 
J1614-2230 specifically in mind, along the line that 
has been developed in anticipation of what is to come 
from the laboratories in construction for the purpose of 
probing cold baryonic matter at high density such as 
FAIR/ G SI and J-PARC. The question we are raising is 
as follows: Is it possible to reconcile the existence of a 
2 Mq neutron star with an EoS that encompasses nor- 
mal nuclear matter, strange hadronic matter, and quark 
matter in a continuous phase diagram? In addressing 
this question, it is, of course, not in our capability to 
make a fully consistent treatment backed by rigorous ar- 
guments. As precisely set forth in the next section, what 
is empirically known in hadron dynamics up to the nor- 
mal nuclear-matter density po that enters into the EoS 
is interpreted in the simplest possible way and extrapo- 
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lated in the given scheme to the unknown higher density 
regime relevant to compact stars, guided by whatever 
available terrestrial as well as astronomical observables. 
What we do is, while sticking as much as is feasible to 
extreme simplicity, develop a chain of arguments which 
are not obviously wrong and could perhaps be quickly 
falsified in the laboratories. In this way, we bypass most 
of the complications invoking "reasonableness." Certain 
caveats inherent in such an approach are pointed out in 
the Appendix. Needless to say, what is badly needed is a 
unified framework in which the different phases involved 
in the process in the compact-star structure are consis- 
tently connected. 

We should mention that there have been a large num- 
ber of treatments in the literature that attempt to take 
into account a variety of relevant degrees of freedom, 
some of which overlap with what we consider here and 
their interactions, as reviewed in Ref. 0. What may dif- 
ferentiate our approach from others is in the - perhaps 
too extreme - simplification of the complexity involved 
in the processes while preserving a certain level of con- 
sistency and that most - if not all - of the assumptions 
we make could be falsifiable. 

Our starting point is the obvious fact that the neutron- 
star matter possesses the n-p asymmetry in its compo- 
sition, the effect of which on EoS is manifested in the 
"symmetry energy," denoted hereafter E sym . x It also de- 
termines the threshold for new degrees of freedom that 
can emerge from weak interactions. The larger the sym- 
metry energy, the more energy the asymmetry between 
neutron and proton composition costs. By introducing 
new degrees of freedom into the system, one can then 
lower the energy of the system and relax the asymmetry 
between neutron and proton. For instance, the energy of 
a pure neutron star matter can be lowered by introduc- 
ing electrons, allowing protons to appear to compensate 
the charge, which, in turn, relaxes the cost of the asym- 
metry owing to the symmetry energy. So unless there 
is any mechanism that keeps the matter from relaxing 
the neutron-proton asymmetry present, the system will 
evolve to a matter in a more stable state by bringing 
in new degrees of freedom. The change of neutron and 
proton fractions of the star matter would require the 
isospin to be violated, which is, of course, not allowed 
by the strong interactions. However, this can happen 
via the weak interactions, where the flavor (isospin or 
strangeness) can be changed. We assume that the evolu- 
tion of the system via weak interactions reaches an equi- 
librium configuration as a ground state of the matter. 
Therefore, the EoS of the star matter in weak equilib- 
rium must be strongly dependent on the symmetry en- 



For those who are not in the field, it should be more aptly called 
"asymmetry energy" because it is the energy that the asymmetry 
between the neutron and proton numbers costs. We follow the 
convention of the community of this field by calling it symmetry 
energy. 



ergy, especially for the emergence of new degrees of free- 
dom accompanied by the change of the neutron-proton 
fraction. One of the most interesting scenarios is that 
kaon condensation can occur inside the neutron star as 
a specific example of a new degree of freedom. We focus 
on this in this paper although other strangeness degrees 
of freedom - when treated differently as specified below 
- can also be important. 

In this work, we investigate the role of kaon conden- 
sation as the principal agent for strangeness in evolving 
to the core of compact stars that exists in the form of 
strange quark matter. Put more precisely, kaon conden- 
sation will play a role of a "doorway" to quark matter in 
the sense precisely stated in the next section. Focusing 
on the symmetry energy, we consider compact stars that 
consist of nuclear matter (NM) in the outer-layer region 
and (1) kaon condensed nuclear matter (KNM) and/or 
(2) strange quark matter (SQM) in the inner region. 
We explore how the structure with KNM and SQM in 
the core region affects the microscopic observables [e.g., 
density-dependent electron (or kaon) chemical potential] 
and macroscopic observables (e.g., mass and radius) of 
the compact stars. 

The outline of this paper is as follows. First we state 
as precisely as possible the strategy that we take for our 
approach. This is the content of Sec. II. The global 
features of kaon condensation and SQM that appear as 
the degrees of freedom extraneous to nucleons is given 
in Sec. III. In Sec. IV, the highly simplified model for 
kaon condensation that we use is described. We also 
discuss the microscopic properties of nucleonic matter 
with kaon condensation and the EoS. The stellar struc- 
ture of this system obtained by integrating the Tolman- 
Oppenheimer-Volkov (TOV) equation is given in Sec. V. 
In Sec. VI, with the given SQM EoS, we study the stel- 
lar structure with SQM driven by kaon condensation. 
The effects of QCD perturbative corrections are also dis- 
cussed. Throughout this work, we set H = c = 1 and 
assume zero temperature. Finally, the results are sum- 
marized and discussed in Sec. VII. 



II. THE STRATEGY 

The system we are interested in involves three basic 
interactions: gravity, weak interactions and strong inter- 
actions. We assume the first two are known, so our focus 
is on the third. We are assuming that solving the TOV 
equation with a given EoS takes care of the gravity in- 
teraction. Both weak interactions and strong interactions 
figure in the EoS. In considering the strong interactions, 
we are concerned - among others - with three possible 
phases, namely, NM, strange hadronic matter, or, more 
specifically, KNM and SQM. Neutron stars, of course, 
have matters other than strong interactions, such as the 
crust which must be considered on the same footing. We 
do not deal with these matters in this work. 

The compact star we are particularly interested in con- 
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sists of three layers: the outer layer of NM, the interme- 
diary layer of kaon condensation, and the inner core of 
SQM. Because we do not know how to formulate the 
problem within the framework of QCD, let us imagine 
having at our disposal a theory that has a hadronic La- 
grangian Chadron applicable up to some density p c , say, 
the chiral transition density or deconfmement density, 
and then above that density, an effective Lagrangian 
C-quark that contains QCD variables, that is, quarks and 
gluons. Apart from the asymptotic density most likely 
inaccessible by compact stars, we know very little of what 
£quark encodes. It could be a color-flavor locked super- 
conducting state, Overhauser state, dyonic crystal, etc. 
We simply assume, as is done by other workers, that it 
can be captured by the MIT bag model with or without 
perturbtive QCD corrections. 

That leaves our principal problem in this paper: 

^hadrow 

In the hadronic sector below p c , several density scales 
must intervene. The first is the NM density denoted po 
and the second is the density at which strange degrees 
of freedom start figuring, in our case, kaon condensation 
at a density denoted p t . There is also - in a specific 
model to be specified below - a possible topological phase 
transition at a density denoted pi/2 lying above po and 
below or above pt- This phase is not understood well at 
the moment, so it is not taken into account in the present 
work. Its possible ramifications in compact-star physics 
are discussed in the Appendix. 

Here is the strategy we take in addressing these differ- 
ent scales. Up to the NM density, the effective hadronic 
Lagrangian Cnn that encodes nuclear physics is mostly 
phenomenological, with the physics involved relying on 
available experimental data. In fact, most of the EoS 
used in the literature have been determined to be con- 
sistent with the available experimental data up to the 
normal nuclear density po, - so all of them more or less 
agree with each other. However, above po, quite arbi- 
trary extrapolations are involved. They are somewhat 
constrained by certain experiments that can probe den- 
sities slightly higher than po- However, the error band of 
the available experimental data is quite large, making it 
difficult to pin down the structure. The presently avail- 
able density regime goes up to ~ 4p with rather large er- 
ror bands. 2 We will pick one particular phenomenological 
parametrization that approximately covers the range of 
the experimental error band up to the density "probed" 
~ 4po- As to the density regime beyond it, we have no 
guidance or control, so we will simply extrapolate it in 
the form we pick. Several potentially important effects 
that are predicted by a given field theoretic framework 
based on hidden local symmetry but cannot be simply 
implemented in the strategy we are taking are mentioned 



2 For an extensive review on the matter on which we heavily rely 
in our work, see [f|. 



in the Appendix. 

At least up to po, the strange degrees of freedom play 
no role, but they must figure at some density above po 
in the guise of kaons and hyperons. In addressing this 
issue, we will take a simplified effective field theory ap- 
proach. Instead of treating both kaons and hyperons at 
the same time, which is quite involved numerically, we 
may integrate out either the kaons or the hyperons in the 
given effective Lagrangian. Now if we integrate out the 
kaons (as well as vector mesons), then the resulting effec- 
tive Lagrangian will consist of a sum of three-flavor [i.e., 
SU(3)f] local four- Fermi (baryon) interactions. Integrat- 
ing out the kaons whose mass becomes small as density 
increases is, in a sense, analogous to integrating out the 
pions from chiral effective field theories for describing nu- 
clear processes whose energy scale is much smaller than 
the pion mass. 4 Doing mean field with this four-Fermi ef- 
fective Lagrangian is equivalent to doing mean field with 
the original baryon SU(3)f chiral Lagrangian in a way 
similar to the Walecka mean-field theory that is obtained 
from four- Fermi chiral Lagrangian (Toj |. 

Alternatively, we can integrate out the hyperons from 
the Lagrangian. This will give us an effective La- 
grangian consisting of kaons coupled to nucleons with 
the couplings determined by the integrating-out proce- 
dure. When treated at mean field (i.e., at tree order), this 
would correspond to the mean-field treatment of chiral 
Lagrangian first put forward by Kaplan and Nelson 
with, however, the coefficients treated accordingly. As 
is elaborated in Sec. IIV A[ there is a hint from lattice 
QCD calculation, that owing to large cancellation be- 
tween higher-order terms, the tree-order treatment sur- 
vives more or less unscathed. Assuming this mechanism 
would render focusing on the kaon sector much simpler 
than looking at the hyperon sector. In this paper, we 
adopt this option. 

To be fair, we should point out that there has been 
a long controversy, still unresolved, on the relevance of 
kaon condensation vs. hyperon condensation in compact- 
star structure. We have no firm argument to support but 
we deem it plausible that both approaches mentioned 
above capture essentially the same physics - and not a 
clear-cut alternative - vis-a-vis to the role of strangeness 
in the dense strange-matter sector of the star. In a more 
complete formulation where both degrees of freedom are 
treated explicitly and on the same footing, there will be 
sharing of effects between the two and most of the ob- 
servables may not be able to disentangle one from the 
other. Future observations of neutrino emission from the 
core might give a hint to the roles of these degrees of 
freedom. 

In sum, our key question is as follows: Within the 



3 This is obviated by the lattice QCD indication that the 
strangeness content of the nucleon is negligible. 

4 This eminently "non-standard" approach to low-energy nuclear 
physics is called "pionless effective field theory" [9|. 
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uncertainties present both in the obesrvables and in the 
parameters of the model, is it possible to develop a sce- 
nario of stellar matter consisting of NM, kaon condensed 
matter, and quark matter with smooth transitions from 
one to the other? 



III. KAON CONDENSATION AND STRANGE 
QUARK MATTER 

In this section, we present a simple description of how 
kaon condensed matter can be linked to SQM. 

In stellar matter, the chemical equilibrium in weak in- 
teraction leads to a condition for the chemical potentials 
of neutron, proton, electron, and muon: 

Pn - Pp = Me = Pu = M- (1) 

The density dependence of neutron-proton chemical po- 
tential difference is determined by the symmetry energy 
in the following form: 

p n - p p = 4(1 - 2^)S(p) . (2) 

Here p is the sum of proton density and neutron density, 
p = pp + p n , and S(p) is the symmetry-energy factor that 
is specified later. 

A. Kaon condensation 

As is discussed below, the key ingredient in our treat- 
ment is the decrease of the effective mass denoted as m* K 5 
of the negatively charged kaon K~ as density increases. 
The m* K is basically a function of mx, p n , Pp owing to 
the kaon-nucleon interactions: 

m* K = u(m K ,p n ,p p ,...). (3) 

The density at which a neutron can decay into a proton 
and K~ via the weak process, n — > p + K~ , 

Pn - Mp = m K, (4) 

determines the threshold of kaon condensation p t . Above 
the kaon condensation, where m* K can be identified as the 
kaon chemical potential px , the chemical equilibrium is 
reached as 

Pn ~ Pp = Pe = Pu = PK = P- (5) 

where 

Pn-p P = 4(l-2^)S(p) + e{K)F(K,p), (6) 



where K stands for the kaon amplitude of kaon condensed 
state, i.e., (K), and F(K,p) is a nontrivial function that 
depends on the neutron-proton chemical potential dif- 
ference which, in turn, depends on kaon-nucleon inter- 
actions. It is a highly model-dependent quantity. The 
charge neutrality condition gives 

Pp = pe + Pfi + &(K)p K . (7) 

Equations ([S]), (O and © are the basic equations to be 
solved to calculate the EoS of KNM. 



B. Strange quark matter 

Strange quark matter can appear as a result of 
confincmcnt-dcconfincmcnt phase transition. This phase 
transition is basically a strong interaction process. Sup- 
pose NM with no strangeness changes over to a decon- 
fined quark phase. There is no room for strange quarks 
in this case because it is the strong interaction that is 
involved. However, in dense compact-star matter, the 
system evolves in weak equilibrium with strange quarks 
in the quark matter through weak interactions. The na- 
ture of interaction for the phase change is, however, dif- 
ferent if it is from a phase in which kaons are present 
to start with induced by the weak interaction. Now if 
kaon condensation has taken place in the system, then 
there is already non- vanishing strangeness in the NM up 
to the phase boundary. In this case, we can imagine the 
confincmcnt-dcconfincmcnt phase transition taking place 
constrained by the weak equilibrium in the kaon con- 
densed matter, leading to SQM. In this subsection, we 
develop the scenario in which kaon condensation joins 
smoothly a SQM. 

We first consider NM (without kaon) and deconfincd 
quark matter (QM). Suppose there is a phase boundary 
between NM and QM at the density p c (for NM) and p® 
(for QM), 6 respectively, in the interior of a neutron star. 
The chemical equilibrium reads 

p n = 2p d + p Ul p p = p d + 2p u . (8) 

The symmetry energies in each phase are related to the 
chemical potentials as 

Pn-Pp = 4(1 - 2x N )S N (p N ) (9) 
Pd-Pa = ±(1-2xq)S q (p q ), (10) 

where x N = P p /pn, xq = p u /pQ, and p Q = p u + p d . 
From the chemical equilibrium 

Pn - Pp = Pd- Pu, (11) 



6 We are taking into consideration that there may be a density 
In this paper, m* K is the kaon energy in medium because we are jump between the NM and the QM at the boundary, such that 

dealing with the s-wave kaon. pj? ^ 3p c . 
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using x Q = Pu/pq = (1 + x N )/3, we get 

(1 - 2x N )S N (p c ) = |(1 - 2x N )S Q (p<?) (12) 

from which we obtain a constraint on the symmetry en- 
ergy at p c , 

Sn(Pc) = \s Q {p% (13) 

for x 7v 7^ 1/2. Just for an illustration, suppose we take 
an expression for the symmetry energy factor of the non- 
interacting quark gas, 

S Q = \{^p Q /2)V*[2V*-l), (14) 

and assume that the critical quark number density, p|?/3, 
is equal to the critical baryon number density, p c \ then 
we will get 

S N ( Pc ) = 24.6(^) 1 / 3 MeV . (15) 
Po 

We may obtain a different constraint on the possible form 
of nuclear symmetry energy, if NM undergoes a phase 
transition into a quark phase at p c . However, that will 
yield a highly nontrivial constraint, which seems not to be 
realized 7 . In Eq. (fl2|) . the simplest solution is xn = 1/2 
irrespective of the specific form of symmetry energy. We 
obtain this result from a different consideration. 

Now consider the case where SQM meets KNM. As 
was done in the previous section, we can obtain the kaon 
threshold density, pt, by the profile of the kaon chemical 
potential. From p t up to the critical density for chiral 
restoration p Cl the weak equilibrium condition reads 

M = 4(1 - 2p p /p)S N ( Pc ) + F(K, p, Pc ). (16) 

Now introduce the Weinberg- Tomozawa term for the 
kaon-nucleon interactions which are introduced in Sec- 
tion [IV2] Since F(K,p,p c ) can be expressed as 

F(K,p,p c ) = pF(K,p c ), (17) 

we have 

pi = 4(1 - 2p p /p)S N ( Pc ) + pF(K, Pc ). (18) 

As the kaon chemical potential - equivalently effective 
mass - p approaches at the critical density, the solution 
x = = 1/2 appears naturally at the phase boundary, 

provided F(K, p c ) =/= 1, which is assumed to be valid for 
the range of density we are concerned with. 



7 In a recent paper by G. Pagliara and J. Schaffner-Bielich Il2( , it 
is argued that the Sq can be three times larger than that given 
by the free quark model if the quark phase is in 2SC, which is, 
of course, a highly correlated system, unlike free quarks. 



The chemical equilibrium (via confinement- 
deconfmcment) reads 

Pn - P P = Pd - p Ul (19) 
Pk- = Ps - Pu (20) 

at the phase boundary. Note that the strange quark is 
required at the boundary, which implies that there should 
be a SQM for p > p c . Because p(= pk) = , we have 
from Eqs. Q3J) and 

p u = p d and ^ s = p u , (21) 

which gives 

Pu — Pd — Ps- (22) 

This is the chemical potential relation for the SQM in 
the masselss limit. In this simple picture, the KNM leads 
naturally to a SQM. 

IV. EFFECTIVE THEORY FOR KAON 
CONDENSATION 

A. The effective Lagrangian 

Following the strategy spelled out in Section [Til we 
consider kaon condensation in the simplest form of an ef- 
fective chiral Lagrangian. Projecting onto the K~ chan- 
nel ignoring all other (pseudo-) Goldstone bosons and 
the nucleons with hyperons integrated out, we write the 
relevant effective Lagrangian in the form [13j 

-C = £-kn + Cnn (23) 

where 

Ckn = d^K-d^K+ -m 2 K K + K- 
+ ^Y> KN {rJn+ptp)K + K- 

+ ■^(n^n + 2 P ^ P )(K + doK- -K-d Q K+) (24) 

£-nn — nhdon + p^idop (Vu' • Vn + Vp^ • Vp) 

2m 

- V NN (25) 

Here the fourth term in Eq. is the well-known 

Weinberg- Tomozawa (WT) term which is constrained by 
a low-energy theorem with / identified with the pion de- 
cay constant /„■ in the matter-free space and Yikn is the 
KN sigma term 

£>kn ~ -(fh + m s )(N\uu + ss\N) (26) 

where m — (m u + md)/2 with the subscripts u, d, and 
s standing, respectively, for up-quark, down-quark, and 
strange quark. In chiral perturbation theory, the WT 
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term corresponds to the 0(p) term and the sigma term 
to 0(p 2 ). In principle, there is no reason why one should 
stop at that order in the chiral counting and also to 
quadratic order in the kaon field as is done in Eq. (|2"4"1) . 
In fact there is a large literature that treats the chi- 
ral Lagrangian to high orders in addressing kaon-nuclear 
physics, but we shall getting into the complexity of such 
sophisticated approaches. 

Although we are unable to make rigorous arguments, 
we can offer several reasons why we can take the sim- 
plest form Eq. (|2"4")l as our basis of analysis. From the 
point of view of renormalization group flow, what mat- 
ters is the WT term with the E term being irrelevant. 
This was shown in Ref. 14]. There it was shown that it 
is the WT term coming with attraction in kaon-nuclear 
interactions that drives kaon condensation. That consid- 
eration does not determine the critical density. It simply 
indicates that kaon condensation is bound to take place 
at increasing density. It also indicates that the A(1405) 
which plays a crucial role for the threshold KN interac- 
tions has little role in kaon condensation transition itself. 
The discussion on these matters - still a highly contro- 
versial issue - is relegated to Appendix. 

That the most important interaction for kaon-nuclear 
physics in the chiral Lagrangian is the WT term is ex- 
ploited in the approach of Ref. [Hj], where all higher 
chiral order terms are ignored and only the WT term is 
inserted as a driving term in a coupled-channel - with 7rE 
channel - equation. In a similar spirit, if one approaches 
kaon condensation from the vector-manifestation fixed 
point [l6j rather than from the matter-free space, one can 
ignore the E term because the quark condensate vanishes 
and retain only the light vector-meson-exchange terms 
(i.e., the WT term) [TtJ ■ The question would arise if 
the E term were as large as thought in chiral perturba- 
tion approach. In the early days, it was thought that the 
strangeness content of the nucleon related to m s {N\ss\N) 
was substantial, say, m s (N\ss\N) > 300 MeV. Thus, in 
ChPT, the E term played the dominant mechanism for 
the condensation. Nuclear density played a role of eat- 
ing into the kaon mass; thus, the "large" E term figur- 
ing importantly in the nucleon structure has an essen- 
tial role. Viewed from the matter-free space, it is the 
E term that is responsible for the relatively low kaon 
condensation density px ~ 3po- This feature is being 
drastically changed by the recent lattice results that find 
that m s (N\ss\N) < 59 MeV which is about one- 
sixth of the old value used in the literature [||. This 
would give a sigma term T,kn % 200 MeV. Furthermore, 
as suggested when viewed from the vector manifestation 
fixed point of view, the sigma term may be strongly sup- 
pressed by density near the condensation point because 
the in- nucleon condensate (N\qq\N) will be suppressed, 
as argued in Ref. [l7|. It seems that the importance of 
the sigma term in the old treatment of kaon condensation 
was overstated or even unfounded. Our point of view is 
that the sigma term Eq. (f2l))) is still quite uncertain at 
high density. 



Next we address the issue of limiting kaon fields to 
quadratic order in Eq. ([24]) for the EoS we are inter- 
ested in. Within the framework, one might attempt to 
calculate a few next orders of chiral perturbation - and 
this has been done in the literature -, but there is no 
compelling reason why a few next orders give reliable 
results. As it appears to be generally the case in ap- 
plying chiral perturbation theory to nuclear and dense 
matter, it may be that going to a few higher-order terms 
that are calculable starting from matter-free vacuum does 
not capture the essential aspect of the process. Indeed 
a recent pioneering lattice QCD calculation [l9] indi- 
cates that this may be the case. In this reference, lat- 
tice QCD simulation of kaon condensation on dense kaon 
systems containing up to 12 AT~'s showed that the prop- 
erties of the condensate are remarkably well reproduced 
by leading-order chiral perturbation theory. While direct 
baryonic background, difficult to implement in lattice cal- 
culations, is still needed, what this surprising result indi- 
cates, rather persuasively, is that there can be substan- 
tial cancelations in higher chiral-order terms that cannot 
be captured by a few perturbatively calculable terms. 
While awaiting for further more quantitative guidance 
from model- independent theory and/or experiments, we 
take the following philosophy in handling the Lagrangian 
Eq. (|2"4")l . Adopting it for our approach - which is essen- 
tially a mean- field one, we simply take the WT as is in 
medium. This is because that term is constrained by low- 
energy theorems. Given the wide uncertainty with the 
sigma term, we will consider T,kn as a free parameter 
constrained by the low value given by the lattice results 
< 200 MeV: It can be as low as ~ 130 MeV depending on 
the 7r-N sigma term as it seems to be required for certain 
nuclear EoS and as high as ~ 400 MeV, as was taken in 
the old days. 

B. The effective Hamiltonian 

Along the strategy described in Sec. |TT] for the kaon- 
nucleon interaction as well as for the nucleonic interac- 
tion, we write the Hamiltonian in the form 

H = Hkn + H nn , (27) 

where 

Ukn = d K-d°K + + [m 2 K - ^-Z KN ]K+K-, (28) 

n NN = Ie° f (£) 2/3 p + v(p) + P (1-2^) 2 s(p). 

(29) 

Here V(p) is the potential energy of NM as a function 
of density p and S{p) is the symmetry energy factor as 
a function of nuclear density. Both S(p) and V(p) are 
model dependent. There are numerous forms for them in 
the literature, all fine-tuned to fit available experimental 
data @. What matters for our purpose is that they be 
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consistent with experiments up to the density that has 
been probed. The data are not precise enough above NM 
density and hence there is a range of parameters that 
are consistent with the error bands given by the data. 
What we will do is to pick one convenient parametriza- 
tion called MDI ( "momentum-dependent interaction" ) 
from the review |8| , which we refer to as "LCK" . It has 
the form 8 



F( V )^ + [18.6-F( V )} 

Po \PoJ 



G(V) 



(2 



2/3 



2/3 



V( P ) 



[m N + ^E° F (P- 
5 \Po 



2/3 



a p 



P 



P 



1 + 7 \Po 



(30) 



(31) 



where a = -298.25MeV, (3 = 244.99MeV, and 7 = 1.21. 
The remaining parameters, F(rj) and G(r]), for this model 
are summarized in Table 1. We take two values for 77, 
i.e., i] = 0,-1. These two values are compatible with 
the band of empirical constraints up to density ~ 4po, 
with 7] = — 1 representing stiffer EoS with the classic 
AKR (Akmal-Pandharipande-Ravenhal) EoS (2(| sand- 
wiched by the two rj values. 9 We see that while consis- 
tent within the error band with experimental constraints, 
they can give different results at higher density in the EoS 
when the strangeness and quark degrees of freedom are 
involved. 



Model 


F(rj) 


G(ij) 




3.673 


1.569 


77 = 


129.981 


1.059 



TABLE I: The parameters for 77 = — 1 and r\ — 0. 
F(rf) is given in MeV unit and G(jf) is dimensionless 
quantity. 



C. S-wave kaon condensation 

Using the chemical potentials defined by 

dU 



(M 



dpi 



(32) 



8 The notation for the variable x in Eqs.(5.1)-(5.5) of Q is changed 
to r) in this paper, because x is used for the proton number 
fraction in the text. 

9 We do not consider r] = 1, which gives what is called "supersoft" 
symmetry energy, which by itself is not consistent with experi- 
mental constraints. See [f|. 




P [P 1 

FIG. 1: The density dependence of the chemical potential 
for LCK with 77 = — 1. The solid, dashed, and dotted lines 
correspond, respectively, to T,kn = 200,300,400 MeV. 



the chemical potentials of proton and neutron are given 
by 



Pp Pp 



p 




2/2 


P 


A + 




p 


P I 



K z 



K 



where 



<-M°=4(l-2^)S(p). 

P 



(33) 
(34) 

(35) 



The amplitude of kaon condensation, K, and kaon chem- 
ical potential, fix , are defined by the ansatz 



Ke 



(36) 



The kaon condensation condition for K ^ is obtained 
by extremizing the classical action, 



2 2 Pn 
m K- P = P— 



2p p 



2/2 



P 



^KN, 



(37) 



which can be solved to get fi or equivalently m* K in 
Eq. ©. 

As stated above, there is only one parameter left unde- 
termined in the kaonic sector. We vary it in what we con- 
sider to be a reasonable range, say, 200MeV < T,kn < 
400 MeV and see how the NM EoS is changed with kaon 
condensation. The density dependence of p, is shown in 
Figs. [Hand [2] for Y, KN = 200,300,400 MeV. In the NM- 
KNM system, the chemical potentials follow these lines. 
The dash-dotted line denotes that of the NM system. 
It can be seen that the chemical potential p begins to 
drop at the kaon threshold density, p t , where the kaon 
amplitude K starts developing non-zero value as shown 
in Figs. [3] and |U We can also see that \i or equivalently 
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FIG. 2: The density dependence of the chemical potential 
for LCK with r\ — 0. The solid, dashed, and dotted lines 
correspond, respectively, to Ekjv = 200,300,400 MeV, as in 

Fig.rj] 




012345678 



P IP ] 

FIG. 3: The kaon amplitude for the LCK model with r\ = — 1. 
See Fig. [T] for the different lines. 




P [P ] 



FIG. 4: The kaon amplitude for the LCK model with r\ — 0. 
See Fig. [2] for the different lines. 



800 




e [MeV / fm 3 ] 

FIG. 5: The EoS for the NM (dashed line) and KNM (solid 
line) with r\ = -1 and Ekat = 130MeV. 



the in-medium effective kaon mass m* K vanishes at some 
high density, which we refer to as critical density denoted 
p c , which could be identified with the chiral restoration 
density. 

The contributions to the energy density and the pres- 
sure from kaon condensation are given by 



2 P v 



K 



Pk = -(m 2 K -p 2 )K 2 



(38) 
(39) 



One can see that the kaon condensation gives a negative 
contribution to the total pressure for px < fn>K- The 



energy density and the pressure of the system are given 
by 



e = V(p)+p(l-2^) 2 S{p) + e lepton 
P 



P = p 



(40) 



2 dV{p)/p + p 2^_ 2 Pp^ ds (p) 



dp 
®{K)P Kl 



dp 



Pi 



epton 



(41) 



where x — p p j 'p. 

To have a qualitative idea of what's going on, we pick 
Y, KN - 200 MeV for r) = and Y, KN - 130 MeV for 
r/ = — 1. We do this because it is the minimum value 
at which the slope of the pressure vs. energy density re- 
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FIG. 6: The EoS for the NM (dashed line) and KNM (solid 
line) with 77 = and Ekjv = 200MeV. 



Model 


Ejfiv 


Pt 




p{pt) 


I/ = -l 


130 


3.74 


13.20 


281.9 


Tj = 


200 


4.27 


8.59 


220.8 



TABLE II: The kaon threshold and critical densi- 
ties for r\ — — 1 and rj = and the Sony's used for 
the TOV calculations (for details see next section). 
pt and p c are given in po unit, and Tikn and A*(p0 
are given in MeV. 



mains positive, so that we can avoid resorting to Maxwell 
or Gibbs construction. We see that the results so ob- 
tained arc not physically reasonable; a higher Y*kn seems 
to be required to be compatible with nature. We list the 
kaon threshold density pt, critical density p c , and chemi- 
cal potential at the kaon threshold pi(pt) for given Tikn's 
in Table 2. 

In Figs. [5] and \5\ are plotted the EoS for Y^kn ~ 
130MeV and - 200MeV. We see that the EoS is soft- 
ened, as expected, beyond the kaon threshold density, 
Pt- 



D. Strange Quark Matter (SQM) Driven by Kaon 
Condensation 



ties 10 , 

Pu = p d = Ps = Pq (42) 
Then the charge neutrality 

2 11 

gP« - gPrf - ^p s = (43) 

is automatically satisfied and there is no need for addi- 
tional leptons. Kaon condensed nuclear matter will nat- 
urally go over to the SQM in SU(3) symmetric phase in 
the massless limit. The EoS of SQM is then given by 

e S QM = a^fj^ + B = 4.83a 4 /9 4 / 3 + B, (44) 
Psqm = 04^5/4 -B = 1.61a 4/ o 4 / 3 — B, (45) 

where B is the bag constant [2l|. Here a 4 denotes the per- 
turbative QCD correction [Q, [22|], which takes the value 
a 4 < 1. The equality holds for SQM without QCD cor- 
rections. In the example of the previous section with 
£kjv= 200 MeV (for 77 = 0), without perturbative QCD 
corrections (i.e., a 4 = 1), the p c comes out to be too high 
so that the bag constant that satisfies the phase bound- 
ary matching condition 

Pknm(Pc) =Psqm(p?) (46) 

becomes negative. Thus, for smaller values of T,kn, the 
kaon-driven scenario may not work. 

If we instead take T,kn — 259 MeV, the pressure 
matching can be satisfied. And we get the critical den- 
sity p c = 6.37po both for 77 = and —1. Using chemical 
equilibrium conditions we get the critical quark number 
densities, p® = 19.26po and 19.11po respectively. One 
can see that pj?/3 are not much different from p c , as 
shown in Table IIIII Given the critical densities, we can 
find a possible set of parameters, a 4 and B 1 / 4 . For ex- 
ample, with a 4 = 0.62, we get B 1 / 4 = 94.71MeV and 
97.25MeV for 77 = and 77 = — 1, respectively. 

The different choice of bag constant with different pa- 
rameters for symmetry-energy, 77, and QCD corrections, 
a 4 , is attributable to the different EoS depending on 77 
and a 4 which results in the different pressures Pknm(Pc) 
for the boundary matching condition, Eq. (|46p . In Figs. [7] 
and [51 we plot the pressure profiles and tabulate the den- 
sities, p t , pt', and p c in Table ITO1 Here p' t {^ p t ) is defined 
by 

PNu{Pt) = Pknm{pv)- (47) 



We now turn to the possibility that kaon condensa- 
tion can drive the dense system to a SQM at the critical 
density p c defined by the condition, pk — 0. For this, 
we assume, for simplicity, massless u, d, and s quarks. 
Eq. (|22p implies that they have the same number densi- 



In a more realistic calculation we have to take into account of 
strange quark mass and electrons and/or muons: p u = 7^ p 3 . 
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FIG. 7: The phase diagram of NM-KNM-SQM for LCK with 
77 = — 1. The solid, dased, and dotted lines are corresponding 
to £14 = 1.0, 0.62, and 0.59, respectively. 
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FIG. 8: The phase diagram of NM-KNM-SQM for LCK with 
rj = 0. The solid, dased, and dotted lines are corresponding 
to 04 = 1.0, 0.62, and 0.59, respectively. 



STELLAR STRUCTURE WITH KAON 
CONDENSATION 



We now turn to a stellar object that has an EoS with 
strangeness. Here we discuss the case where kaon is con- 
densed but without SQM and the case where kaon con- 
densed and with strange quark core. 



Model 


Pt 


Pf 


Pc 


P?/3 


t/= -1 


3.05 


4.95 


6.37 


6.37 


rj = 


3.76 


5.03 


6.37 


6.42 



TABLE III: The characteristic densities of NM- 
KNM-SQM with the LCK models for r\ = -1 and 
77 = with Ekjv — 259MeV. The values in the last 
column are the quark critical densities which are 
equivalent to the baryon number densities, respec- 
tively. 
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FIG. 9: The schematic phase diagram for NM-KNM. 



A. Double-layered (NM+KNM) structure 

First we consider a simpler case, where neither 
Maxwell nor Gibbs construction is needed at the kaon 
condensation phase transition as shown in a schematic P 
vs p curve (Fig. [5]) . In the framework of our model with 
the LCK parametrization for EoS, the maximum values 
for ^kn turn out to be - 130 MeV and ~ 200 MeV for 
i] = — 1 and respectively. For those values of T,k n > we 
cannot have pressure balance for SQM unless we allow a 
negative bag constant. Hence up to the critical densities 
in these parameters, one finds NM and KNM from the 
outer layer to the core part. We use Eqs. P0")) and (|4Tj) 
with Eqs. (J3S]) and ([39]) to integrate the TOV equation 



dM 
dr 
dP 
dr 



Aner 2 , 



GMe 



P 

1 + - 

e 



47rr 3 P\ / 2GM\ 



where M(r) is the mass enclosed inside the radius r. 

The results are given in Figs. 10 and 11 They illustrate 
how softening the EoS with kaon condensation affects the 
mass and the size: It lowers the mass and increases the 
size from the case of NM at the same central density. 
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density 



FIG. 10: M-R sequence for LCK with 77 = -1 and T, KN = 

130 MeV. The solid and dashed lines are for the cases with FIG. 12: The schematic phase diagram for NM-KNM-SQM. 
NM-KNM and NM, respectively. 
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FIG. 11: M-R sequence for LCK with 77 = and E K n = 
200 MeV. The solid and dashed lines are for the cases with 
NM-KNM and NM, respectively. 



shall simply take the maximum mass in this case to be 
given by the mass with the critical density, p c , of Table 
lllll taken as the central density of the star. However, with 
the stiffer symmetry energy with r\ = — 1, the maximum 
mass is obtained by the gravitational instability condi- 
tion at the central density much smaller than the critical 
density but not far from the kaon threshold density. This 
feature can be seen in Table IIVI 



Model 






r .Vf. -V.,,,.,,.,. 


NM (77 = -1) 


2.23 


11.41 


6.1 


NM+KNM (77 = -1) 


2.10 


12.33 


5.2 


NM (77 = 0) 


1.90 


9.81 


8.1 


NM+KNM(77 = 0) 


1.74 


9.61 


8.5 



TABLE IV: The calculated structure parameters 
for various systems (see text for details). The units 
are Mq, km, and po for mass, size, and central 
density, respectively. 



However, because we use the relatively small values for 
XifAr, the effect of kaon condensation is not so significant 
(e.g., about ~ 0.2M Q reduction in the maximum mass). 

In Table IIVI are summarized the stellar structure pa- 
rameters of our results for both the NM and the NM- 
KNM systems. It is interesting to note that the two 
parametrizations for the symmetry energy, with 77 = 
and -1, give rather different features of M vs. central 
density relations. For the softer symmetry energy, 77 = 0, 
the mass M increases with the central density beyond 
the kaon condensation threshold density going up to the 
critical density. It implies that the star is gravitation- 
ally stable up to the critical density. Beyond the critical 
density, we have no models to calculate the EoS. So we 



B. Triple-layered (NM+KNM+SQM) structure 

If we take a larger value for Hkn, say, 260 MeV, it 
is possible to consider a triple-layered structure shown 
schematically in Fig. Q2] consisting of NM, KNM, and 
SQM from the outer layer to the core part. In contrast 
to the case considered above, there will be an unstable 
region between p t and pt> , for which we have to resort to 
Maxwell or Gibbs construction. In the Maxwell construc- 
tion, it turns out to be impossible to match the chemical 
potential, which signals instability. In the Gibbs con- 
struction, however, there appears a mixed phase of NM 
and KNM. In this work, we take a rather simple ap- 
proach, namely, allow the discontinuity of density and 
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FIG. 13: The M-R sequences for LCK with 77 = -1. The 
dotted line denotes NM. The dashed-dotted line between A 
and B denotes the double-layered (NM-KNM) system. The 
solid lines, a, b and c, denote the triple-layered (NM-KNM- 
SQM) system with the QCD corrections with 0,4 = 0.59, 0.62 
and 1, respectively. 



FIG. 14: The M-R sequences for LCK with 77 = 0. The dotted 
and dashed-dotted lines are the same as in Fig. 1131 The solid 
lines, a, b and c, denote the triple-layered (NM-KNM-SQM) 
system with the QCD corrections with 0,4 = 0.59, 0.62 and 1, 
respectively. 



chemical potential by assuming that NM with the den- 
sity pt changes into KNM with the density p' t at the phase 
boundary denned by 

P{p t )=P{p v ). (49) 

The resulting mass-radius relations that follow from the 
TOV equation are plotted in Figs. [T3] and [TU As ex- 
pected, kaon condensation with higher Scat's would lead 
to smaller masses. However, when the central density 
becomes higher, an SQM driven by kaons appears at the 
core. No sharp change due the emergence of SQM is 
observed, which implies that the kaon-driven SQM tran- 
sition is a rather smooth transition in this scenario. 



Model 


04 




-^1 M = M maX 


P {c) \ 


r? = -l 


0.62 


1.99 


10.60 


12.3 


T) = 


0.62 


1.85 


9.86 


13.5 



TABLE V: The maximum masses, M max , radii, 
R, and central densities, pj c ) for the three-layered 
structure with the QCD corrections, 0,4, = 0.62 for 
both of 77 = —1 and respectively. The units of 
the stellar structure parameters are the same as in 
Table [TV] 

The results given in Table IVl highlight the importance 
of perturbative QCD corrections in viable three-layered 
structure. We find that with the EoS of SQM without 
perturbative corrections, i.e., a 4 = 1, the results (the 
solid line c in Figs. [13] and H~4| are similar to the previ- 
ous case with the soft symmetry energy (with 77 = 0) of 



the kaon condensed matter, which is present up to the 
critical density. Here the mass increases with the central 
density even beyond the kaon threshold density to reach 
the maximum mass, which is slightly greater than the 
mass at the kaon threshold. However its detailed struc- 
ture is quite different from the latter with Y^kn = 200 
MeV, where there is no room for the SQM. If, however, 
we include the QCD corrections, the maximum mass lies 
beyond the critical density, and the core is in the form 
of SQM. With the stiffer symmetry energy with 77 = — 1, 
the maximum mass of ~ 1.99Af Q can be obtained by 
the gravitational instability condition at the central den- 
sity p = 12.3 po. With the softer symmetry energy (with 
r\ = 0), the maximum mass of ~ 1.85M0 can be obtained 
at the central density p = 13.5 po, as shown in solid lines 
b in Figs. 13 and 14 and in Table E] with a 4 = 0.62. 
For different set of parameters with 04 = 0.59, we get 
slightly larger maximum masses, 2.03Af Q and 1.90Af Q , 
respectively, for 77 = —1 and 77 = 0, as shown with solid 
lines a in Figs. 13 and 14. This is a rather high density 
and a more sophisticated treatment of SQM than the sim- 
ple MIT bag model with perturbative corrections (e.g., 
color superconductivity) may be called for and change 
the maximum mass. 



VI. SUMMARY AND DISCUSSION 

We have proposed a scenario in which dense compact- 
star matter is driven smoothly to an SQM by kaon con- 
densation at the density at which the kaon chemical po- 
tential [Lk — rn* K becomes negligibly small. It is based 
on the main assumption of the paper that kaon conden- 
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FIG. 15: Pressure vs. density (p) for NM-KNM (solid line) 
and NM-KNM-SQM (long-dashed line) with LCK (r) = -1) 
and considered parameters in previous sections compared 
with the analysis by Ozel et al. (shaded area)[23||. The dot- 
ted line denotes NM, which is given for comparison with NM- 
KNM and NM-KNM-SQM. 



sation does occur at a finite density and that the kaon 
effective mass should become very small at a higher den- 
sity, which, however, should not be larger than that of 
compact star core. We show that the compact star with 
an NM-KNM-SQM structure is possible with the param- 
eters that are not excluded by theory or phenomenology. 
We admit that our calculation is highly model-dependent 
without the benefit of constraints by QCD or by reli- 
able models, but it is simple and could be easily fal- 
sified by experiments. All three states of matter, that 
is, NM, kaon condensed matter, and QM, are treated on 
the same level of accuracy. For kaon-nuclear interactions, 
we take the simplest chiral effective Lagrangian at tree 
order, which seems to be supported by the available in- 
medium lattice QCD with the Y,kn term being taken as 
the only unknown quantity. As for nuclear interactions, 
we take LCK's empirically parameterized form for the 
energy density and symmetry energy. There are poten- 
tially serious effects that are missing in the treatments 
and that could drastically change the scenario. They are 
discussed in Appendix. 

We find that the maximum mass of compact stars with- 
out kaon condensation with the LCK parameters can be 
twice the solar mass, which is consistent with the recent 
observation. With kaon condensation, the EOS becomes 
soft and the maximum mass gets reduced, as expected. 
With E K jv = 200 MeV (and 77 = 0), an SQM driven by 
kaon condensation in our scenario is impossible because 
the pressure matching at the SQM boundary gives a neg- 
ative bag constant. However with the somewhat larger 
^kn = 260 MeV, an SQM driven by kaon condensation 
can be formed with the bag constants B 1 / 4 ~ 188 and 
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FIG. 16: The same as Fig. [T5] for LCK (77 = 0). 
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179 MeV, respectively, for 77 = — 1 and 7/ 
maximum mass can be as large as 2Mq for 77 = 
A in Fig. Q21 where only NM is relevant. For 77 = 0, a 
triple-layered structure is possible, but without pertur- 
bative QCD corrections in SQM, the maximum mass is 
far below 2Mq, as given by the solid line c in Fig. [14] 
When suitable QCD corrections are included in the EoS 
of SQM, smaller bag constants are required, 12 97.25 and 
94.71 MeV, respectively, for 77 = —1 (04 = 0.62) and 
?7 = (04 = 0.62). It is also found that the transition is 
smooth as far as the mass change is concerned. This in- 
dicates that the kaon condensation as we described here 
can be considered a "doorway" to a quark-star matter in 
the core. 

What we learn from our analysis made in this paper 
is that the nature (stiffness or softness) of the EoS in 
the nucleon sector represented in the LCK model by the 
different 77 values (even within the range allowed by ex- 
periments), the sigma term Tikn, the bag constant, and 
QCD corrections is highly correlated. Indeed, for the 
stiffer EoS with 77 = —1 and the smaller Y^kn ~ 130 
MeV, a two- layer structure of NM and KNM can be made 
compatible with the mass vs radius constraint given by 
the 1.97 M Q star. With the larger E KN ~ 260 MeV, 
even the three-layer structure of NM, KNM, and SQM 
with an appropriate bag constant is possible beyond the 
critical density, 6.37po- And perhaps most significantly, 



Note that for the case of the kaon-condensation-driven SQM, 
SjfJV — 259 MeV is fixed for rj's. However, the bag constants 
need to be different to match the pressure condition. 
The bag constants needed for the case with perturbative correc- 
tions may appear to be too low in comparison with the values 
needed in vacuum phenomenology. This may not be a cause for 
worry, however. Because the part of the gluon condensate locked 
to the quark condensate is expected to "melt" across the chiral 
transition point, the smaller bag constant effective in the SQM 
phase appears to be more natural. 
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the maximum mass of a neutron star with the three-layer 
structure can be made consistent with the recent obser- 
vation. There is, of course, nothing new or surprising in 
this sort of correlations in nuclear processes. However, it 
points for a reliable treatment of the M vs. R ratio to the 
necessity of treating all of the variety of different aspects 
of the phenomenon on the same footing. This point can 
be illustrated by the following observation. With a rela- 
tively stiff EoS with nucleons only such as that typified 
by the one constructed in Ref. [20| (which is comparable 
to the LCK's with r\ = —1), the critical density for kaon 
condensations can be pushed to too high a value to be 
of relevance to the process, say, > 7po [24|. However, as 
discussed in the Appendix, an effective field theory con- 
sideration based on hidden local symmetry suggests that 
at a certain high density approaching the density regime 
where quark degrees of freedom emerge, the repulsion 
that pushes up the critical density for kaon condensa- 
tion may get strongly suppressed as one approaches the 
critical density. This means that the NM and KNM sec- 
tors have to be treated in the same framework. Similar 
remarks can be made for the SQM sector. 

In our treatment, strangeness was brought in by kaons. 
It could have been done by hyperons. As argued in Sec- 
tion |lTl the two ways at the mean-field level should cap- 
ture approximately the same physics. They may not be 
representing alternatives. We expect that they will give 
qualitatively same effects on the stellar structure. 

The results we obtained are compared with the analysis 
by Ozel, et al. [23( shown in Figs. [T5] and [TH] as shaded 
areas. In view of the variety of intricate mechanisms left 
out in our treatment, it is not obvious that one can do 
much better with more detailed model calculations. As 
stressed, what is needed is a consistent field theoretic 
approach that can account also for the effects of the type 
discussed in Appendix. We are only at the beginning of 
such an endeavor. 
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cal effects on the kaon propagation in dense matter and 
possible binding of kaons in finite and infinite systems, 
the fate of w-meson-mediated repulsion at high density. 
These issues are potentially relevant to the problem ad- 
dressed in our paper. To make a reliable treatment tak- 
ing all these considerations into account would require 
a fully consistent field theoretic treatment and such a 
treatment is presently lacking in the literature. Most of 
the treatments found in the literature are a piece-meal 
patching of different ingredients with no obvious connec- 
tion between them. We are unable to do much better 
here either. What we can do is just to point out how the 
various ingredients can figure in kaon-nuclear and high 
density physics. In so doing, however, we have in mind a 
specific field theoretic framework based on HLS [r| and 
its generalized form - which by itself is consistent but 
incomplete - although we do not have the formulation 
to allow us to do the calculation explicitly within that 
framework. 



(A) The role of A(1405) 

In considering the possible mechanisms for or against 
the formation of dense kaonic (finite) nuclei, the A(1405) 
has figured prominently. There is a humongous contro- 
versy on whether and how A(1405) enters in generating 
strong or weak A(1405)-nuclear potential. We gloss over 
this controversy here because it does not seem to be es- 
sential to our approach. What we would like to describe 
here is that, regardless of what the A(1405) does near 
the threshold AT~-nucleon interactions and how it gen- 
erates the A(1405)-nuclear potential, it is irrelevant to 
our approach to kaon condensation. If one were to try 
to arrive at the putative kaon condensation by starting 
from elementary kaon-nucleon interactions, the structure 
of the A(1405) would seem to matter. In fact, it has 
been suggested that even if there is a strongly attractive 
kaon-nuclear interaction, kaons do not condense macro- 
scopically [251 ] . There is a difference in perspectives on 
this issue. 

In Sec. I IV A[ in addressing kaon condensation, we 
made no account of the presence of A(1405). This is 
because, as was first pointed out in Ref. [26|], the pre- 
cise location of A(1405) is irrelevant to the phase tran- 
sition involved. This point is depicted schematically in 
Fig. El 13 

There Z? _1 (lu) gives the dispersion formula for the kaon 
frequency w. 14 Given that kaon condensation occurs far 



Appendix 

In this appendix, we address a variety of issues that are 
not taken into account in our work: The role of A(1405) 
in dense kaonic nuclei and kaon condensation, topologi- 



We are grateful to Chang-Hwan Lee for allowing us to use this 
figure. 

We are, of course, assuming that we can treat A(1405) as a quasi- 
particle in medium. If this were not the case, the spectral func- 
tion could be more complicated. Even so, the qualitative feature 
of our argument will still hold as suggested by the RG argument 
of Ref. pj. 
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One has to start 
below A (1405) pole ! 



T> 




without A(1405) 

with A( 1405) 



D- 1 (w) = w a -M&-II(w) 



FIG. 17: (Color online) Irrelevance of A(1405) in kaon con- 
densation. Mf( is the effective kaon mass. 



below the threshold, that is, with M% far below the 
A(1405) pole, the A(1405)-nuclear interaction is irrele- 
vant in the RGE sense of Ref. [HJ so that it makes no 
difference whether A(1405) is there or not. What this 
tells us is that kaon-nuclear interactions in dilute system 
have no influence on the presence of the phase transition, 
though they could influence its location. 



(B) Effect of topology change in dense matter on 
EoS 



describing nuclear and dense matter [27J ■ No such treat- 
ment that is reliable is available to date. One reason why 
it can be taken seriously in the absence of a credible jus- 
tification is that it involves a topology change. As such, 
its qualitative feature could be robust, more or less free 
of model dependence. In fact, this feature presents an 
advantage. Normal NM is most likely a liquid with its 
stability assured by the Fermi-liquid fixed point - and not 
a solid. However, it is possible that the matter is solid at 
high density, as suggested by N c arguments. Now what 
we are observing on the crystal is a change of topology, so 
the properties that are anchored on topological structure 
could remain valid in what we are looking at although 
the density is not high enough. 

That the Skyrmion-to-half-Skrymion transition takes 
place at some high density is more or less model- 
independent within the adopted framework, although de- 
tail features will depend on models. In Ref. (28|, such a 
calculation was performed with the Skyrme Lagrangian 
implemented with the "soft" dilaton that enters in the 
trace anomaly of QCD that is assumed to melt across 
chiral restoration [29J. From this model, certain qualita- 
tively robust features emerge. At pi/ 2 where a Skyrmion 
fractionizes into two half-Skyrmions, the quark conden- 
sate (qq) (in QCD variable) which is given in the Skyrme 
model by TrU (where U is the chiral field) vanishes on 
the average in one unit cell; thus, chiral symmetry ap- 
parently "restored" but with a non-zero pion decay con- 
stant /* 7^ 0. In this phase, quarks are still confined 
in hadrons. Such a putative chiral symmetry restored 
but quark-confined phase has not been observed in con- 
tinuum models. If confirmed to be viable, however, it 
will have an important impact on various phenomena, as 
described below. 



(Bl) Skyrmion-half-Skyrmion transition 

When dense matter is simulated by putting Skyrmions 
on crystal lattice, which is potentially a powerful tool for 
looking at strong-coupling dynamics of baryonic matter 
at high density 27], one discovers that there is a sort of 
phase transition at a certain density above NM density 
that we denote as pin > po from the Skyrmion matter - 
which corresponds to the usual NM - to a half-Skyrmion 
matter populated by half-soliton configurations. In that 
phase, chiral symmetry is putatively restored but quarks 
are still confined. Although the critical density pi/ 2 can- 
not be pinned down, it presumably lies between 1.3 and 
2 times the NM density. This transition is quite novel 
and has never been predicted before and there are no ex- 
perimental data to confirm or refute this prediction. It 
could perhaps be vindicated or falsified in future exper- 
iments at such laboratories as FAIR/GSI, J-PARC, or 
RIB machines in construction. 

There is a strong indication that hidden local sym- 
metric Lagrangian, especially with an infinite tower as 
recently discovered, will be more reliable and more pre- 
dictive than the Skyrme model with the pion field only for 



(B2) Anomalous kaon binding 

One intriguing consequence of the topology change is 
"anomalous" binding of K~ in dense Skyrmion mat- 
ter. As stressed above, a powerful treatment for dense 
baryonic matter could be made with a three-flavor effec- 
tive Lagrangian with hidden local symmetry, but this is 
numerically involved and has not been worked out yet. 
Though less reliable, one can instead take the Skyrme La- 
grangian as an effective one that results when all vector 
mesons are integrated out. Given such a Skyrmion La- 
grangian, one can treat the property of a kaon in dense 
medium in terms of the kaon propagating in the back- 
ground of Skyrmions put on a crystal lattice. This was 
worked out in 28| . The main agent for binding the kaon 
to the Skyrmion in this description is the Wess-Zumino 
term in the SU (3) / Lagrangian. This attraction is the 
analog to the Weinberg- Tomozawa term in the chiral La- 
grangian Eq. (|2~4")l . However, there is an additional mech- 
anism for kaon-nuclear attraction not apparent in the La- 
grangian that comes from the Skyrmion-half-Skyrmion 
phase change. The result of Ref. [28[ is summarized in 



16 



Fig. [TBI where the effective kaon mass m* K is piotted as a 
function of density p/po (denoted in the figure as n/no) 
for two different values of the (free-space) dilaton mass. 15 
What we note is that the location of pi/ 2 is remarkably 
insensitive to the dilaton mass, but m* K varies strongly 
with density at p\/ 2 and beyond. Thus the following 
warning: What is given in Fig. [TS] should be taken with 
caution. 




FIG. 18: m *K± vs - p/po in dense skyrmion matter which 
consists of three phases: (a) (qq) / and f£ 7^ 0, (b) (qq) — 
and f^^=0 and (c) (qq) = and /* = 0. The parameters are 
fixed at \/2e/^ = m p — 780 MeV and dilaton mass m x = 600 
MeV (upper panel) and m x = 720 MeV (lower pannel). 

There are a several noteworthy observations to make 
on this result. One is that there is a sudden - and siz- 
able - decrease in m* K at p\/ 2 . Such an "anomalous" 
drop cannot take place without topology change. Thus, it 



15 The precise value of the dilaton mass - as well as the microscopic 
(QCD) structure of the dilaton — is unknown. This is the part of 
the puzzle on the nature of scalar mesons in low-energy hadron 
physics. The values used here are what are considered to be 
reasonable values for the mass of the scalar that enters into mean 
field theory of NM. 



cannot be present in standard many-body or field theory 
treatments of nuclear and dense matter in which topol- 
ogy does not figure. What we learn from this observation 
is that if one were to produce a compact kaonic-nuclear 
system that can reach the density corresponding to pi/ 2 
- which should be above, but not too high above, the 
normal matter density, there the system would undergo 
a drop in kaon mass that is not present in the standard 
treatment. 16 It is perhaps significant that a similar "dis- 
continuity" in nuclear tensor forces takes place at the 
same transition density p 1 / 2 mentioned below. 

The other interesting observation is that the kaon mass 
drops down to a small value, with the kaon disappear- 
ing, at the point where /* falls to zero, presumably sig- 
naling the deconfined phase. In our scenario of kaon- 
condensation - quark-matter continuity discussed in the 
text above, we are imagining an analogous process with- 
out, however, the apparent discontinuity of the magni- 
tude that depends on the dilaton mass that is seen in 

Fig. mi 



(B3) Tensor forces, "new BR scaling" and E sym 

Another potentially important consequence of the 
topology change described above is on the nuclear tensor 
forces. It will therefore affect the EoS of compact-star 
matter, particularly on the symmetry energy E sym 
This aspect, highly relevant to our main discussion, is, 
however, not explicitly taken into account in the paper. 
It will be done in a future presentation 

In the effective field theory approach to nuclear forces, 
in which the relevant degrees of freedom are nucleons, 
pions and vector mesons, the nuclear tensor forces are 
highly sensitive to in-medium properties of the p meson, 
specifically its mass and its coupling to nucleons. The 
p tensor force comes with the sign opposite to that of 
the pion tensor, so their contributions to the tensor force 
tend to cancel as density increases because the strength 
of the p tensor force increases at increasing density. This 
is because the parameters of the p meson scale according 
to what is known as BR scaling 32], where the p mass 
drops at increasing density following roughly the quark 
condensate in medium. There are evidences for this scal- 
ing in a variety of nuclear processes, among which the 
most striking one is the quenching of the Gamow- Teller 
matrix element responsible for the long life-time of 14 C, 
i.e., the famous 14 C dating [33]. There may be alterna- 
tive explanations for the quenching of the Gamow- Teller 
matrix clement involving many-body potentials or many- 
body correlations. However, within the given effective 
field framework used in Ref. [33[, it provides a strong 
support that the BR scaling does hold up to NM density. 



A similar anomalous drop is proposed in the work of Akaishi et 
al [30| . They referred to it as "enhancement effect." 
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However, there is no evidence either for or against the 
BR scaling of Ref. [32j above the NM density which is 
pertinent for neutron stars. 

In [3l| . it is shown that the Skyrmion-half-Skyrmion 
transition can drastically modify the BR scaling [32| at 
the topology-change density pi/2- Skipping details, we 
simply state what was found in that reference. A sub- 
tle scaling change occurs at p l / 2 owing to the topology 
change which makes the p tensor force whose magni- 
tude increases up to the density px/2 turn over and start 
decreasing, ultimately get suppressed more or less com- 
pletely at a density near two to three times that of NM. 
This topology change turns out to leave the pion tensor 
unaffected. Therefore at some high density above P1/2, 
only the pion tensor remains. In the model calculation of 
Ref. [3l|, the complete suppression of the p tensor takes 
place at about 3poi in the vicinity of kaon condensation. 
In terms of the symmetry energy, a cusp appears at pi/2- 
Whether this is real or an artifact of the crystal structure 
is not clear. In any event, this phenomenon can have a 
dramatic effect on processes in baryonic matter at high 
density, in particular on the symmetry energy in the den- 
sity regime that has not been accessed by the presently 
available experiments. Future experiments will provide a 
check for this prediction. 



(B4) Dilaton limit and suppression of the 
short-distance repulsion 

Consider a baryon HLS Lagrangian implemented with 
a dilaton scalar introduced via the trace anomaly of 
QCD. We call it dHLS. Such a theory at mean field is 
known to be able to describe NM fairly well provided 
the parameters of the HLS Lagrangian are endowed with 
suitable BR scaling pyj. The question raised is as follows: 
How can one use this Lagrangian to describe the chiral 
phase transition that is to take place at high density, say, 
p c ? This question was answered in (35j . 

While it is not known how to approach systemati- 
cally as function of density the chiral restoration point 
with this Lagrangian, it was suggested by Beane and van 
Kolck I36| how to approach Weinberg's mended symme- 
try limit [33] by taking what is called the "dilaton limit." 
In that limit, the dilaton-implemented HLS Lagrangian 
goes over to Gell-Mann-Levy linear sigma model which 
can then be used to describe the chiral phase transition in 
the standard way. This limiting process involves a trans- 
formation of a chiral singlet - which the dilaton is - to 
the fourth component of chiral four-vector, namely, er of 
the sigma model. Apart from the symmetries involved, 
this may be caused by a sort of level crossing in terms 
of quark configurations. In the quark-gluon language of 
QCD, the dilaton could be a complicated combination 
of four-quark q 2 q 2 wavefunction and glueball wavefunc- 
tion that is low-lying at low density while the a is the 
qq component which lies higher. The level crossing may 
occur such that at high density near the chiral restora- 



tion, it is the latter that comes down low in energy. The 
key point here is that as density approaches that of chi- 
ral restoration, the dilaton limit transforms dHLS to the 
linear sigma model endowed with the mended symmetry. 
It was found in Ref. p35| that the dilaton limit in dHLS 
corresponds to taking, among others, 

9a = gv -> 1 (Bl) 

where gA is the axial coupling constant in the neutron 
beta decay and gv is related to the V-nucleon coupling 

as 

gvNN = ffws(l - 9v) (B2) 

where V = (p, w) are the lowest-lying vector mesons and 
ghis is the hidden gauge coupling. 

There are a variety of phenomena that could follow as 
consequences of the approach to the dilaton limit and 
influence strongly high-density physics. What is closely 
relevant for this discussion is that as one approaches the 
dilaton limit (|B1|) . the V-nucleon coupling decreases and 
goes to zero at the dilaton limit. Owing to the vector 
manifestation of HLS [101, the hidden gauge coupling 
going proportionally to the quark condensate (qq) near 
chiral restoration will also go to zero. However the vec- 
tor meson mass also goes to zero as ~ ghis, so the ratio 
gvNN/ m v ~ constant. This means that as density in- 
creases into the regime where the dilaton limit is valid, 
interactions that involve vector-meson exchanges will be- 
come weak independently of the approach to the VM 
fixed point. In particular, the NN repulsion at short 
distance attributed to u> exchanges will be suppressed 
proportionally to (1 — gv) 2 - This would make a dras- 
tic change to the conventional thinking of the hard core 
repulsion that has been understood up to date to be op- 
erative at high density. Furthermore if one were to as- 
sociate short-distance repulsion involving three nucleons 
with u) exchanges, the suppression factor would also go 
like (1 — gv) 2 and so be equally suppressed. From what 
density this suppression becomes active is not known, 
but if this mechanism is confirmed to be viable, it will 
certainly affect the presently available EoS' at densities 
beyond what has been "measured." Among various phe- 
nomena, as noted in the Discussion section, such a sup- 
pression mechanism could invalidate the argument of Ref. 
24] - based on short-range repulsion in nuclear interac- 
tions - against low critical density for kaon condensation. 
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